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1/Q^ corrections in DIS fragmentation functions 
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We present results of recent calculations [ [jj, H of power corrections to single-hadron inclusive momentum 
distributions in the current hemisphere of the Breit frame in DIS. Though results are presented for both singlet 
and non-singlet scattering contributions, we emphasise the fact that the former is dominant in the kinematical 
region studied at HERA. 



1. INTRODUCTION 

The work described in the following article is 
largely motivated by the recent experimental in- 
terest in the study of hadronic final state mo- 
mentum distributions in the current hemisphere 
of the Breit frame of DIS. The observable in ques- 
tion is similar to the hadronic energy distribution 
in e'^e~ annihilation and is related in an analo- 
gous manner to the universal parton to hadron 
fragmentation functions. However in order to 
compare theoretical predictions for these DIS mo- 
mentum distributions accurately with the HERA 
data, we must take into account the existence of 
1/Q^ non-pcrturbativc effects (power corrections) 
from both the singlet and non-singlet processes. 

We wish to study the distribution F'^{z), de- 
fined by 



F\z;x,Q^) 



d'^cr'^ 



dxdQ^dz I dxdQ"^ 



(1) 



for a given hadron species /i, as a function of the 
variable z — 2ph • q/q^, where ph is the four- 
momentum of the resultant hadron and q the 
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photon four-momentum. Since the fragmentation 
products of the remnant of the nucleon are ex- 
pected to be travelling in directions close to that 
of the incoming nucleon, i.e. in the 'remnant 
hemisphere' pz < 0, we consider only hadrons 
produced in the 'current hemisphere' Pz > 0. 
Such hadrons are expected to be fragmentation 
products of the scattered parton. Thus z takes 
values < z < 1. The other quantities x, Q" 
have their usual DIS meanings. 

One has to consider power corrections to both 
the numerator and the denominator of Eq. (|^) 
each of which have a leading power correction 
that goes as l/Q^. The denominator is just 
the familiar DIS cross section which can be de- 
composed in terms of transverse and longitudinal 
structure functions for which the results are pre- 
sented in Refs. [|j and [^. For the numerator one 
has an identical decomposition in terms of gener- 
alised {z dependent) structure functions to which 
we also compute power corrections. We present 
here, results for both singlet and non-singlet con- 
tributions to F^[z;x,Q'^). The calculations de- 
scribed here were done using the dispersive tech- 
nique which we briefly comment on below. 



2. CALCULATIONS AND RESULTS 

The calculation for the non-singlet piece fol- 
lows the usual dispersive method first described 
in [ pj . In this method power corrections corre- 
spond to non-analytic pieces in a fake gluon mass 
e, in the 0{as) calculation of a given observable, 
represented by the characteristic function J-'{e). 
Importantly, the normalisation of the result is in 
terms of hopefully universal parameters, which 
turn out to be various moments of a universal 
strong coupling as{fi'^), defined for all scales ^^. 

In order to calculate power-behaved contribu- 
tions to QCD observables involving more than a 
single renormalon chain, like the singlet sector 
here, we need to generalise the standard disper- 
sive treatment, the details of which can be found 
in references [ 0, H. Here we just mention that in 
cases involving for instance two gluons, like the 
forward scattering amplitude of the DIS singlet 
graphs, we need to introduce two dispersive vari- 
ables (gluon masses) and our characteristic func- 
tion J-'{ei, 62) is now an 0{a^) quantity. However 
it turns out that the answer can be expressed in 
terms of functions J- which are obtained from the 
characteristic function by setting one of the gluon 
masses to zero while retaining the other. Hence it 
is still the non-analytic terms in a single variable 
e which correspond to power corrections. 

In the present case (DIS fragmentation func- 
tions) we find the following types of non- 
analyticity which translate into 1/Q^ corrections 
as below: 



T ^ ct log e 



and 



T r^ ic'elog^e 



&F 
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jDi 



log- 



ics 



Q2 "Q 



where Di and D2 are defined by 

/•OO 



(2) 



(3) 



(4) 



1 1"°° 
\0gD2 = — / d^l'' log ^l^ [2as5as - {5a,f] .(5) 
^i Jo 



and 



represents a non-perturbative modification to the 
standard perturbative form of thr coupling. We 
emphasise the fact that in the above formula a^"^ 
stands for the perturbative expansion of a^ip'^) 
in terms of as{Q^), where this expansion is trun- 
cated at whatever order the perturbative result 
for a given observable has been computed. Hence 
in most cases of current interest this would be at 

We calculate the singlet and non-singlet contri- 
butions to the fragmentation functions, retaining 
terms up to 0{e) non-analytic at e = 0. We find 
terms that diverge as log e and log e arising from 
collinear splitting; the gluon mass here behaves 
as a regulator, and the divergence is factored into 
the scale dependence of the parton distributions. 

The 1/Q^ corrections arise from terms propor- 
tional to e log e and e log e. In the singlet case we 
obtain the correction 



h Di TjiCp 



SF" = 



g2 (27r)2 






(6) 



where Kt{z, x) and Kl{z, x) are shown in figure 
g and y is the usual DIS variable. Note that if 
Di is positive, we would expect to see a large 
negative power correction at small z. 

In non-singlet scattering we find the relative 
correction 



5F" 



A 



Fh Q 



7^H{z,x) 



(7) 



Sa^ifj, ) = as{^ ) - a^'^(^ ) 



where H{z, x) is shown by figure ||. This gives 
a large positive correction at large z, but only a 
small correction at small z. 

The detailed calculations and algebraic expres- 
sions for results can be found in [ |l|, |j . 

3. CONCLUSIONS AND PROSPECTS 

The results presented in this article clearly in- 
dicate that at small values of the fragmentation 
variable z, one would expect the singlet contribu- 
tion to be dominant. Moreover at small values of 
the Bjorken variable x, we would anyhow expect 
the singlet scattering mechanism to start becom- 
ing important due to the rise in the gluon density. 

The parameter Di that appears in the singlet 
result has to be phcnomcnologically determined 
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Figure 1. Values of H{z,x) for non-singlet frag- 
mentation. Dashed, dot-dashed and solid curves 
are quark, gluon and total fragmentation; the two 
sets of curves are for x = 0.3 (upper) and 0.1 
(lower) . 




Figure 2. Graphs showing Kt{z) and Kl{z) 
singlet fragmentation at fixed values of x. 



(extracted once and for all from the data). Re- 
sults are relatively insensitive to the other param- 
eter D2. Unfortunately however, before one can 
compare the predictions made here with the ex- 
perimental data one needs to clarify several is- 
sues. 

From a technical viewpoint one has to learn 
how to treat the singlet piece as a gluonic con- 
tribution rather than an 0(0^) non-singlet con- 
tribution. In other words, the result presented 
here was obtained by treating the incoming gluon 
(that scatters via qq production off the photon) as 
an internal line radiated off an incoming quark. 



This is the only way to obtain the power cor- 
rection. However to reinterpret this as a singlet 
contribution we have to find a meaningful way to 
factorise off the quark to gluon splitting and then 
convolute with the gluon density rather than the 
quark density. 

Then there is also the question of doing im- 
proved perturbative calculations. At present the 
NLO program CYCLOPS [ § is used to obtain 
the perturbative prediction. However in actual 
fact, at small values of z we would expect terms 
that vary as log(l/z) to spoil the convergence of 
the perturbative series. Hence the need arises 
for resummed perturbative estimates which take 
these terms into account at all orders. 

Another relatively minor consideration is that 
for the present results we calculate the distribu- 
tion differential in z, the hadron longitudinal mo- 
mentum fraction, rather than the energy fraction, 
since it turns out to be complicated to analyti- 
cally compute the latter in DIS. Hence when com- 
parisons are made with data one must make sure 
that we are using the same fragmentation variable 
although qualitative features can be expected to 
be similar for both choices. 

In fact at this conference we saw experimen- 
tal results for the energy distributions where the 
deviation from NLO perturbation theory was in- 
deed large and negative as we expect at small z 
[ . Remarkably the data can be extremely well 
described using a simple power correction ansatz 
of Dokshitzer and Webber also presented here. 
However the ansatz is not a theoretical predic- 
tion in any strict sense but rather an educated 
guess at the shape of the correction. Hence for 
a full theoretical picture we need to understand 
the problems mentioned above. Once this is done 
we can proceed with extracting Di from the data 
and using it in fits to other singlet observables. 
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